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Abstract. A generic 'chirp' of the form h(i) = A(t) cos <j>(t) can be closely 
approximated by a connected set of multiscale chirplets with quadratically- 
evolving phase. The problem of finding the best approximation to a given signal 
using chirplets can be reduced to that of finding the path of minimum cost in 
a weighted, directed graph, and can be solved in polynomial time via dynamic 
programming. For a signal embedded in noise we apply constraints on the path 
length to obtain a statistic for detection of chirping signals in coloured noise. In 
this paper we present some results from using this test to detect binary black hole 
coalescences in simulated LIGO noise. 



1. Introduction 

Despite having achieved unprecedented sensitivities, experiments for laser interfero- 
metric detection of gravitational waves such as LIGO pQ face significant challenges, 
not least of which is the problem of detecting unmodelled or poorly modelled sources 
of gravitational waves. For detecting the inspiral of a binary system, the standard 
technique is matched filtering using a bank of templates parametrised by the compo- 
nent masses of the system. For low-mass binaries, the time evolution of the inspiral is 
well-modelled by post-Newtonian approximations, however for high-mass binaries the 
models are considerably less certain [2]. Furthermore, as the binary mass increases, 
the spin of the two bodies becomes a significant factor in the evolution of the signal [3] . 
A complete description of a binary system including the spin of both bodies requires 
17 parameters, making the set of templates to be searched over infeasibly large. Even 
when some parameters are neglected, estimates of the number of templates needed to 
detect, for example, spinning extreme mass ratio inspirals using a space-based detec- 
tor such as LISA range from 10 15 -10 40 templates [H[5]. Methods have been proposed 
to reduce the number of templates required, such as by using detection template fam- 
ilies which cover the expected range of gravitational wave signals [H [7] , but these still 
require ~ 10 5 templates [S]. 

Template methods for detecting binary coalescence events mostly focus on the 
inspiral component or the ringdown component [9\ and do not attempt to match 
the merger component, believed to be a major contribution to the gravitational 
signature for black hole coalescences. Modelling the inspiral and ringdown is relatively 
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straightforward, whereas modelling the merger requires robust techniques for solving 
the full Einstein equations numerically under extreme conditions. Much progress has 
been made in achieving this goal but the problem is far from solved [TU] E] ■ 

A number of potential gravitational wave signals are of short duration (less than 1 
second) and are collected under the heading of burst sources. These include events such 
as supernovae, the final stages of binary black hole coalescence, and other potential 
sources of gravitational waves such as gamma-ray bursts. Generally, models for these 
sources are either non-existent or insufficient for constructing matched filters, and we 
must rely on non-parametric methods. Various methods for detecting bursts have 
been proposed [12]-[l7], and some have been applied to interferometer data [15] . 

In this paper we apply a non-parametric detection scheme called the best path 
(BP) test introduced in [19] to the detection of binary black hole coalescences in 
simulated LIGO noise. The terminology comes from the study of weighted graphs and 
refers to the path between two vertices of a graph which is of maximum total weight, 
subject to a constraint on it's length. The BP test is applicable to the detection of 
quasi-periodic signals of the form 

h(t) = A(f) cos 0(f) (1) 

where the amplitude A(t) varies slowly with time and the unknown phase <p(t) obeys 
some regularity conditions. Signals of this form have a well-defined instantaneous 
frequency f(t) — 4>(t)/2it (to avoid confusion, we note that there is an unrelated 
method called the Fast Chirp Transform which is applicable to the detection of signals 
of the form |T]) where the phase function is known 20 ). 



2. Chirplet path pursuit 

Given detector output 

u{t) = n(t) + ph{t) (2) 

where n(t) is Gaussian coloured noise with 2-sided power spectral density S(f), we 
seek a test statistic which will discriminate between the two hypotheses 
H : p = 

H x :p^Q. (3) 

The null hypothesis is that the data is pure noise, while the alternative is that the 
data contains a chirp- like signal of the form JT]), normalised with respect to the inner 
product derived from S(f), 

f ii*(f)v(f) ... 

{u > v)= L^ur df - (4) 

The parameter p may be interpreted as the expectation of the SNR, 

SNR=^-. (5) 

rms (n, n) 

Note that a 1-sided PSD is more commonly used in the literature, equivalent to 2S(\f\). 
We use the 2-sided PSD here to simplify the discretised form of 

Locally, chirps with smoothly-varying phase have a very simple structure. Over 
short times their frequency evolution is approximately linear. For longer duration, 
local approximations can be joined together so that the instantaneous frequency of 
the signal is approximated by a piecewise linear function. In the following we outline 
the methodology for obtaining a test statistic via chirplet path pursuit - details may 
be found in [T9l. 
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2.1. Multiscale chirplets 

Consider a signal on the interval I = [0,T). The preceding discussion suggests we 
should examine functions which will correlate well locally with signals of the form (JTJ) . 
Our detection method uses a dictionary of normalised multiscale chirplets of the form 

c SJ , a , & (t)cxe l2 ^ t+fct2 / 2 ), teI s jQI (6) 

that is, a collection of chirplets supported on intervals I s j and parametrised by length 
scale s, location j, initial frequency a and chirp rate b. The intervals are taken to 
be dyadic of the form I s j = [j2~ s T, (j + 1)2~ S T]. Here s = 0, 1,2, . . . represents a 
scale index and defines the length of the dyadic interval. The dictionary has elements 
of various durations, locations, initial frequencies and chirp rates. It is convenient to 
think of a chirplet as a line segment a + bt supported on I s j in the time-frequency 
plane. 

Our test statistic is constructed by looking for a connected 'path' of chirplets in 
the time-frequency plane that gives a good overall correlation with the signal. To 
achieve this we notionally discretise the time-frequency plane and consider points 
(U, fk) as vertices in a directed graph. The frequency intervals may be chosen as 
convenient - for example, to coincide with bins of a discrete Fourier transform. Fixing 
a time- frequency discretisation also fixes the the discretisation of the chirp parameter, 
since we think of chirplets as arcs connecting vertices of the graph supported on dyadic 
intervals. Using the FFT we can quickly calculate the local correlations \(u,c s j : a : b)\ 2 
of u(t) with elements of the chirplet dictionary, which we use as the weights of the 
arcs connecting each vertex in the graph. Given a connected, non-overlapping chirplet 
path P = {ci, C2, . • • , c p } supported on a partition V = {Ii,l2, ■ ■ ■ , I p } of I the total 
weight of the path is | (u, c p ) | 2 . A description of our discretisation scheme may be 
found in the Appendix. 

Simply maximising ^2 p \{u,c p )\ 2 ov er all chirplet paths will naively overfit the 
data. In the limit of small chirplets, such a statistic would simply fit u(t) rather than 
a hidden signal. Instead we use a multivariate statistic obtained as the solution of the 
optimisation problem 

T £ * = max ^ |(u,c p )| 2 subject to |P| < I. (7) 

v 

Here £ is a constraint on the path length ie. the number of chirplets in the path. To 
be adaptive, we calculate T £ * for several different path lengths, I E L = {£i,£2, ■ ■ ■}■ 
While there are a vast number of possible paths, using a variant of Dijkstra's algorithm, 
calculating T/ reduces to a constrained dynamic programming problem which can be 
solved in 0(\L\ x #arcs) [21]. The number of arcs depends on such things as choice 
of discrete frequencies and chirp rates, but is typically not more than N 2 log 2 N. 

Since T£ is a multivariate statistic we use a multiple comparison rule for rejecting 
the null hypothesis 22]. Given data u(t) we test Hq at false alarm probability a using 
the following procedure: 

1. For each £ € L, calculate T£ and find the corresponding p- value under Hq, pf. 

2. Compare the minimum p- value p* — miiiipi with the distribution of minimum 
p- values under H . 

3. If p* is small enough to lie in the a-quantile of the distribution, reject Hq - we 
conclude a signal is present. 
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In this procedure we are choosing the ordinate of the multivariate test statistic that 
gives the greatest evidence against the null hypothesis. We then compare this p-value 
with what one would expect under the null hypothesis. Although there do not exist 
analytic expressions for the distributions of Tj* and the minimum p-value, we can 
estimate them using Monte Carlo simulations. We call the best path (BP) statistic. 
As an example, Figure [1] shows the path obtained for an inspiral signal in white noise. 



Time-frequency image of 




Figure 1. Best path found for a binary inspiral signal with total mass 16 Mq in 
white noise, indicated by the dashed curve. Vertical lines delimit the support of 
individual chirplets in the path. Notice that the BP test uses long chirplets when 
the frequency is changing slowly, and short chirplets when it is changing rapidly. 



3. Simulations 



3.1. Noise model 

To estimate the statistical power of the BP test we have studied the detection of certain 
gravitational wave signals in simulated LIGO noise. Discretely sampled Gaussian noise 
is produced via the following method. We generate two sequences of white noise a&, 
then construct a discrete Fourier representation of an instance of coloured noise hk 
using the PSD as follows: 
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(8) 



By construction, the inverse DFT n% is real Gaussian noise with PSD Sk- The PSD 
used is the polynomial fit given in [23l Table 5]. This fit is only valid for frequencies 
above the LIGO-I seismic wall frequency f s = 40 Hz. Seismic noise renders region 
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below f s inaccessible to gravitational wave searches. For the purposes of simulation, 
we mimic high-pass filtered data by rolling off S(f) below 20 Hz. When calculating 
the BP statistics we only search over paths with instantaneous frequencies above f s . 
At this time we have not included non-Gaussian features such as instrumental bursts 
in our noise model. 



3.2. Signal model 

Since the object of the exercise is to detect 'real' gravitational waves, we will use as 
our test signals a collection of physically realistic waveforms for binary black hole 
coalescence. We use a modification of the method in 12 to model a complete 
coalescence waveform. The signal consists of an inspiral component, a merger 
component, and a ringdown component. While the inspiral and ringdown models 
are reasonable, the simulated merger should not be taken to be physically realistic. 
Instead, it is meant to approximate the overall time and frequency characteristics of 
a real merger. 

The test signals are parametrised by the total mass M = mi + m-i of the two 
bodies and the symmetric mass ratio r/ = m 1 m2/M 2 . The full waveform is obtained 
by combining the components in such a way that the instantaneous frequency and 
amplitude are continuous up to first derivatives: 

( A ins P(i)cos0 ins P(t) i<0 
h(t)= I A mc ^ c (t) cos </) mcisc (t) 0<t<t m (9) 
[ A Tin &(t) cos ct) Tin s(t) t m <t. 

Here we have arranged for the inspiral component to end at t = and the merger 
component to end at t — t rn . Following |24j we take the merger duration to be 
t m = 50M/M Q x T Q . 

For the inspiral component of the signal we use the non-spinning 2PN 
approximation for the phase in the form given by [251 eqn. 15.24]. For amplitude 
we use the leading order (ie. Newtonian) expression given in [25], eqn. 15.27-28]. For 
simplicity we average over orientation ((,, (3) and sky position to obtain 

2/3 

(10) 



A insp {t) = ^. T&cllM 



ttTqM f lnsp (t) 



M, 



5 D M Q 

where D is the distance to the source. 

We model the inspiral component from the time the instantaneous frequency 
enters the sensitive band of the detector above f s up to the commencement of the 
merger component. Deciding where the boundary between inspiral and merger lies is 
somewhat arbitrary. We follow [24] in making the transition at the point where post- 
Newtonian approximations begin to break down. It is convenient to fix this transition 
at t = 0. A conservative estimate [24] is that errors in the 2PN approximation become 
significant when the instantaneous frequency reaches 

fo = ^M X 4100 HZ (U) 
so we set the coalescence time t c of the inspiral in J35J ec L n - 15.24] by solving 
/ insp (0)=/o- 

The ringdown component is assumed to be an exponentially damped sinusoid with 
constant frequency / rm s as given in [551 ec l n - 18.3]. Our amplitude model, adapted 



Gravitational wave detection using multiscale chirplets 



6 



from [25], is 



A rins {t) 



A Tqc M 



e -Kf"*{t-t m )/Q 



(12) 



where a is the dimensionless spin parameter, Q — 2(1 — a) 45 is the quality factor, 



and e is the fraction of M radiated as gravitational waves during the ringdown. The 
factor of 1/V207T in (|12jl comes from averaging over orientations and sky positions. 
This is essentially the same amplitude model as given in [25j eqn. 18.5]. 

Our inspiral component has been arranged to terminate at t = 0, with ringdown 
commencing at t = t m . Since no analytic models exist for the merger component, we 
fit the amplitude and phase functions to bridge the gap between inspiral and ringdown. 
Assuming that the merger waveform is of the form |T]), a simple way to connect the 
inspiral and ringdown waveforms is to require that the amplitude be continuous to first 
derivatives, and the phase to be continuous up to second derivatives (thus ensuring 
that the instantaneous frequency is continuous up to first derivatives). This gives four 
conditions that must be satisfied by / mer s e (£) and A merse (t) at t — and t = t m , so we 
model J mer s e (t) and ^4 mergc (t) by cubic polynomials. Since we also require the phase to 
be continuous at t = 0, we obtain </> merge (i) from the anti-derivative of / mor s°(£) with 
an appropriate constant of integration. We note that phenomenological templates for 
coalescing binaries have recently become available [27], however our waveforms are 
qualitatively very similar, and for testing purposes it is convenient to know the exact 
form of the instantaneous frequency and be able to set the precise times of transition 
from inspiral to merger to ringdown. Phenomenological templates will be examined 
in future work. 

3. 3. Choice of signal parameters 

To test detection efficiency we used signals of length N = 512, N = 1024 and N = 2048 
sampled at 2048 Hz. Signals of roughly this duration are produced by BBH systems 
with total mass in the range 20-50 M G . As most models for the ringdown waveforms 
assume equal mass binaries, we will only consider this case. The masses used were 
mi = m,2 = 22.5, 15 and 10. Motivated by recent numerical experiments [101 [11], we 
take a = 0.7 and e = 0.01. While the procedure for producing a merger waveform is 
crude, it does produce a signal with frequency and amplitude characteristics similar 
to those seen in numerical relativity simulations. Figure [2] shows the strain and 
instantaneous frequency for these binary coalescences at a distance of 1 Mpc for the 
M = 45 and 30 M cases. 

4. Results 

To use the BP test we first need the distribution of Tt under Hq. There is no analytic 
expression for the distribution of the BP statistic so we have used a Monte Carlo 
simulation to estimate them. As our test signals have different lengths we generated 
three null distributions, one for each N. In each case we generated 10 5 instances of 
simulated LIGO noise and calculated T/ for each of them with chirplet path lengths 
£ drawn from the set L = {1, 2, 4, 8, 16}. These random trials give an approximation 



-,1/2 



A = 4 



Q [1-0.63(1 -a) - 3 ] 



(13) 
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Figure 2. h(t) and instantaneous frequency for a binary coalescence with masses 
(a) mi = m,2 = 22.5 Mq and (b) m\ = m2 = 15 Mq at a distance of 1 Mpc. 



to the distributions of under Hq for each t. Using our empirical distributions we 
can estimate the p- value for an observed T£ . 

To test detection efficiency, we first constructed normalised test signals using the 
model described in Section I3T2"! For each p — 7,8,9, 10, 11 and 12 we generated 10 5 
instances of noise and injected the signal at that level. The BP statistic Tjf was 
calculated, as was p*, and we determined the detection probability for a given a 
by counting the number of p* < a. Figure [3] gives the detection probabilities as a 
function of a (the Receiver Operating Characteristic curve) for p — 8, 10 and 12. For 
comparison, we also give an ROC curve obtained using matched filtering to detect the 
signal. For these curves p has been chosen to give a good match to the ROC curve 
obtained via the BP test. From this it can be seen that the BP statistic is about 
half as sensitive as matched filtering. Since the distance D to the source is inversely 
proportional to the overall signal amplitude, we can consider the BP test to have a 
seeing distance about half that of matched filtering. 

In Figure [3] we also give the detection probability as a function of p (or 
equivalently, inverse distance to the source) for a = 0.05, 0.01 and 0.001. The 
corresponding distances at p = 10 are D — 100 Mpc for M = 45 Mq, D = 80 Mpc 
for M = 30 Mq and D = 65 Mpc for M — 20 Mq. This shows that, for example, at 
a false alarm probability of a — 0.001 we can see an event out to ~ 100 Mpc with a 
false dismissal probability of about 10%. Note that since we have averaged the signal 
amplitude over sky positions and orientations, an optimally aligned and positioned 
source could be detected much farther away. 

In the above comparison we are injecting a known signal into noise and using 
the same (normalised) signal as our template for matched filtering. Real signals in 
interferometer data will have unknown parameters, and a bank of templates using 
discrete values of the parameters (mass, spin etc) is needed to cover the range of 
physically plausible coalescences. Since a real signal has parameters drawn from a 
continuum there will usually be some degree of mismatch between the signal and 
templates in the bank. As such, the comparison above is very conservative in 
comparing the BP test with the most favourable matched filter detection scenario, one 
which is unlikely to be attained in practise. A more realistic benchmark is obtained 
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Figure 3. Detection probability as a function of p and false alarm probability 
for BBH coalescences with total mass (a) M = 45 Mq (b) M = 30 Mq and (c) 
M = 20 Mq. In each case we give an estimate for the p which gives a similar 
curve using matched filtering. 



by examining the performance of the BP test when the signal parameters are chosen 
at random from a range of values. Here we present a comparison of the BP test with 
detection via a bank of templates, and with another method employed in searches for 
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unmodelled signals, the excess power statistic [13] . 

We first created a bank of templates using discrete values for the parameters. 
Although our complete signal model contains a large number of free parameters, for 
simplicity we chose to only vary mi, mi and a, as these have the greatest effect on 
the waveform. For the same reason we have used equal spacing in all parameters, 
rather than attempting to construct a template bank spaced to give equal overlap 
between adjacent templates. While methods exist to construct optimally-spaced 
template banks, our templates have the additional complication of including merger 
and ringdown components. 

For each of the signal lengths N = 512, 1024 and 2048 we generated a bank of 
normalised templates using the criteria that 

1. The range of masses mi and m.2 is chosen so that the length of the signals range 
from N/2 to N samples. 

2. The spin ranges from a = 0.18 to 0.98. 

3. The spacing between masses and spins is chosen so that the minimal match of a 
signal with parameters drawn from the range of parameters is at least 0.97. 

For each N we then generated 1000 test signals with mass and spin parameters 
drawn at random from the appropriate range, and injected them into simulated LIGO 
noise with p = 10. The resulting data was used to calculate a BP statistic for each 
segment. Comparing the BP statistics with the empirical null distribution as above, 
we obtained the ROC curves shown in Figure 3J Searching for the same signals via 
matched filtering, we found that the ROC curves matched well when the signals were 
injected with p around 6.5 - in other words, the BP test sees about 2/3 as far as the 
template bank. 

We performed a similar analysis using the excess power statistic, which is optimal 
when the only known features of the signal are the duration and bandwidth [13J. 
The excess power statistic is simply the power (u, u) calculated using Q where the 
integration is performed over the bandwitch of the expected signals, taken to be 40- 
1024 Hz in this instance. Under H this has a x 2 distribution with degrees of freedom 
twice the number of frequency bins. Using the excess power statistic, we found that 
the ROC curves matched those of the BP test in Figure Q] well when p was around 14. 

5. Conclusion 

Chirplet path pursuit has previously been shown to be effective at detecting a broad 
class of chirp-like but otherwise unmodelled signals in coloured noise 19J. In this paper 
we have demonstrated that the method can be successfully applied to the problem of 
detecting test signals with similar characteristics to those expected from binary black 
hole coalescence. The method is able to detect a range of signals of modest strength 
hidden in simulated LIGO noise, and exhibits somewhat better statistical power than 
the excess power test. 

As with other methods for detecting bursts, in real LIGO noise there is the 
difficulty of distinguishing genuine gravitational wave signals from instrumental and 
environmental events. For matched filter searches the \ 2 discriminator can be used 
to reject signals that do not have the correct distribution of power across frequency 
bands, however this requires that the gravitational waveform be known [55]. This 
discriminator is not applicable to chirplet path pursuit since the signal is not known 
and we do not impose any assumptions on the distribution of power. Instead we 
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N = 512, mass range 20.5 -30.4 M Q N = 1024, mass range 13.5-20.4 W Q 




(a) (b) 



N = 2048, mass range 9.5 - 1 3.4 W Q 




(c) 

Figure 4. Detection probability as a function of false alarm probability for 
random signals in the mass ranges (a) 20.5-30.4 Mq (b) 13.5-20.4 Mq and (c) 
9.5—13.4 Mq. In each case we give an estimate for the p which gives a similar curve 
using matched filtering with a bank of templates and the excess power statistic. 

would rely on the methods being employed in current searches: requiring events to be 
coincident across multiple detectors, vetoing events based on environmental channels, 
and testing if waveforms measured in different detectors are consistent [T5] . 

As expected for a non-parametric method, chirplet path pursuit is not as sensitive 
as matched filtering using a template bank, nevertheless our comparison shows that 
the method has similar effectiveness to matched filtering for a signal that is roughly 
1.5 times as strong. Significantly, since the method is sensitive to a wide range of 
chirp-like signals, an exact model of the signals to be detected is not necessary. This 
makes the method particularly of interest in situations where the signal is unmodelled 
or poorly modelled, as is the case for the late inspiral and merger components of 
intermediate mass black hole coalescences. 
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Appendix 

In this appendix we present the scheme used for calculating chirplet coefficients of 
discretised data. The data u[n\ = u(nAt) is discretely sampled at N = 2 s intervals of 
duration At. Notionally, this discretises the time-frequency plane into points (£«,/&) 
where U = iAt and fk — kAf = k/(NAt). Points in the time-frequency plane are 
considered to be vertices in a directed graph where the weight of the arc connecting 
two vertices is given by the local correlation of u(t) with the corresponding chirplet. 

Consider chirplets supported on the interval [0, 2~ S T). At scale s, < s < S this 
interval has length N s = 2~ S N samples. While there are many ways to discretise 
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chirplets on this interval, it is convenient to choose the spacing of the frequency 
parameter to correspond with the bins of a discrete Fourier transform, and choose 
the spacing of the chirp parameter so that at the end of the interval the instantaneous 
frequency has changed by a whole number of bins. Thus our dictionary of chirplets is 
indexed by scale index s, frequency index k and chirp index I, and the (unnormalised) 
discrete chirplet is given by 

c s ,fc,iW = e l2 ** ,[nl 0<n<N s , < k < N/2 (A.l) 

where the phase is 

TL V? 

<j) s k An] = k— + I— — — . (A.2) 

vs.b.u j N 2NN S y ' 

The discretised instantaneous frequency is 

TL 

^ S ,k,i[n] = k + l—. (A.3) 

-< » s 

Such a chirplet has initial frequency kAf Hz and rises to frequency (k + l)Af Hz 
at a rate of lAf /(N s At) Hz s _1 . Since we only deal with real signals, the range of 
the chirp index I is chosen to restrict the chirplets to non-negative frequencies up to 
Nyquist, thus -k < I < N/2 - k. 

In general, the inner product |4]) for a noise process with covariance matrix X is 
where u* is the conjugate transpose of u. For our noise model the Fourier 
matrices F mn = e -^ mn l N diagonalise S, and so (u,v) — u*F*D~ 1 Fv = u*D~ 1; i) 
where D = diag(oo, af , . . . , a%_x) and g\ — (\nk\ 2 ) are the eigenvalues of E. 
Calculating (u, c) for a chirplet supported on a dyadic interval I s j = [j2 _s T, (j + 
1)2~ S T) is equivalent to calculating the inner product of u(t + j2~ s T) with a chirplet 
supported on [0, 2~ S T). As the time index of the first sample in I s j is jN s , let 
u s ,j — (u]jN s ],u[jN s + 1], . . . , u[(j + l)N s — 1]) be the samples of u(t) restricted to 
I s .j. Then to find (u,c) we pad u s j and c to length N with zeroes and FFT. In 
discrete form, the inner product then reduces to 

At ^ u* [n]~c[n] 

(u ' c) — £^r (A - 4) 

n = 1 1 

where S[n] — S(nAf). If c has indices k, I then after normalising, | (u, c)\ 2 /\(c, c)| 2 is 
the weight of the arc connecting (t jNa , f k ) to (ty+i)jv, i fk+l)- 

To calculate the BP statistic we must find the total weight of connected, non- 
overlapping chirplet paths in the time-frequency plane starting at t — and ending at 
t = T . To keep the number of arcs manageable we further restrict our chirplet paths 
to those supported on a recursive dyadic partition (RDP) of / constructed using the 
following definition [3D] : 

1. The trivial partition V = {1} is an RDP. 

2. If V = {I\,l2, ■ ■ ■ ,I p } is an RDP, then so is the partition obtained by splitting 
any interval Ij into two adjacent dyadic intervals. 

This means that, for example, that {[0, 1/4), [1/4, 1/2), [1/2, 1)} is a recursive dyadic 
partition of [0, 1), but {[0, 1/4), [1/4, 1)} is not. The total weight of a chirplet path 
P = {ci, C2, ■ ■ ■ , c p } supported on V = I2, ■ ■ ■ , I P } is then 

\\ c Pt c P/\ 



